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Inflaton or Curvaton? 
Constraints on Bimodal Primordial Spectra from Mixed Perturbations 
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We consider Cosmic Microwave Background constraints on inflation models for which the primordial 
power spectrum is a mixture of perturbations generated by inflaton fluctuations and fluctuations in 
a curvaton field. If future experiments do not detect isocurvature modes or large non-Gaussianity, 
it will not be possible to directly distinguish inflaton and curvaton contributions. We investigate 
whether current and future data can instead constrain the relative contributions of the two sources. 
We model the spectrum with a bimodal form consisting of a sum of two independent power laws, with 
different spectral indices. We quantify the ability of current and upcoming data sets to constrain 
the difference An in spectral indices, and relative fraction / of the subdominant power spectrum 
at a pivot scale of fco = 0.017 Mpc -1 h. Data sets selected are the WMAP 7-year data, alone 
and in conjunction with South Pole Telescope data, and a synthetic data set comparable to the 
upcoming Planck data set. We find that current data show no increase in quality of fit for a 
mixed inflaton/curvaton power spectrum, and a pure power-law spectrum is favored. The ability 
to constrain independent parameters such as the tensor/scalar ratio is not substantially affected by 
the additional parameters in the fit. Planck will be capable of placing significant constraints on the 
parameter space for a bimodal spectrum. 
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I. INTRODUCTION 

Current astrophysical data are consistent with a primor- 
dial spectrum of cosmological density perturbations con- 
sisting of a power-law spectrum, with spectral index close 
to scale invariance. This is consistent with the predic- 
tions of the simplest single-field models of cosmological 
inflation. In the single-field scenario, fluctuations in the 
field responsible for inflation (the inflaton) generate a 
spectrum of curvature perturbations in the early uni- 
verse. In this case, the spectral index of the resulting 
perturbation spectrum is determined by the shape of the 
inflaton potential. However, other mechanisms for gen- 
erating curvature perturbations are possible, for example 
the curvaton scenario PHI] , in which a massless "specta- 
tor" field sources isocurvature perturbations during the 
period of inflationary expansion. The isocurvature per- 
turbations are subsequently converted to curvature per- 
turbations when the curvaton field decays, which is as- 
sumed to happen after the end of inflation and reheating. 
Like single-field inflation, the simplest curvaton models 
predict a power-law spectrum of perturbations, but in the 
curvaton case, the shape of the spectrum is completely 
unrelated to the inflaton potential. 

Even when a curvaton is present, however, inflaton 
fluctuations will still generate curvature perturbations 
[5] . While it is typically assumed that perturbations from 
the inflaton are subdominant to the curvaton spectrum, 
there is no physics requiring this to be so. It is in princi- 
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pie possible for both power spectra to be of similar ampli- 
tude - in fact, suppressing the inflaton spectrum requires 
fine-tuning of the inflaton self coupling even beyond the 
A ~ 10~ 14 already required by Cosmic Microwave Back- 
ground (CMB) normalization, so it is reasonable to ex- 
pect that inflaton-generated perturbations will be signif- 
icant even when a curvaton is present. An important 
distinction between inflaton and curvaton-sourced per- 
turbations is that, while the curvaton three-point corre- 
lator, the so-called bispectrum, can be sizable [Sit?! , that 
of the inflation is suppressed by slow-roll parameters and 
is therefore totally negligible [H [9] . A future detection of 
non-Gaussianity (NG) in the CMB will therefore point to 
a light field like the curvaton as the source of the curva- 
ture perturbation [TO] . Additionally, a residual isocurva- 
ture component in the primordial spectrum will indicate 
the presence of the curvaton [TT] . However, suppose that 
future experiments do not give any sign of NG or isocur- 
vature modes. How can we know if the curvature pertur- 
bations come from the inflaton or if they are created by 
a different mechanism, e.g. by converting isocurvature 
perturbations into adiabatic ones? This is an empirically 
relevant question, since single field models of inflation will 
be constrained under the assumption that the curvature 
perturbations originated from the inflaton. The result is 
that a given set of CMB observables can be mapped to 
a wide range of different inflaton potentials, each with a 
different, unresolved curvaton contribution. This inver- 
sion was studied in Ref. |llj . where it was found that 
the resulting degeneracy in the inflaton potential param- 
eter space is large. While additional observables beyond 
the adiabatic and tensor two-point correlators are ulti- 
mately needed to disentangle the inflaton and curvaton 
contributions, in this work we examine whether current 
and future data can instead shed light on their relative 



contributions to the overall curvature spectrum. 

In this paper we therefore consider CMB constraints 
from current and upcoming data on a mixed infla- 
ton/curvaton primordial power spectrum, which we 
model as a superposition of two uncorrelated power laws 
with independent spectral indices. The paper is struc- 
tured a follows: Sec. [XT] specifies the model studied and 
defines the parameters fit to CMB data. Section |III| 
presents constraints from the WMAP 7-year (WMAP7) 
data, WMAP7 plus South Pole Telescope (SPT) data, 
and simulated Planck-quality data consistent with cur- 
rent constraints. Section HVl summarizes our conclusions. 



II. BIMODAL POWER SPECTRUM 

In this work we consider scalar perturbations formed 
from the superposition of two uncorrelated power law 
spectra, 



P(k) = A x 



A, 



n 2 — 1 



(1) 



Bimodal spectra of this form can arise if both the in- 
flaton and an additional decoupled field contribute to 
the overall curvature spectrum. As a specific physical 
model to motivate Eq. (IT]), we consider mixed inflaton 
and curvaton fluctuations. The curvaton, a, is a weakly 
coupled scalar field that is relatively light during infla- 
tion (m 2 <C H 2 ). It is decoupled from the inflationary 
dynamics and over-damped by the expansion, 3H a ~ 0, 
until m 2 < H 2 . After inflation, the curvaton is free to 
roll to the minimum of its potential, V(a), where it un- 
dergoes oscillations during a radiation dominated phase. 
The fluctuations in the energy densities of the radiation 
and curvaton field generate curvature perturbations, ipi, 
given in terms of the gauge invariant expression, 



G = -ip- H — 
Pi 



(2) 



The variable £.; reduces to the curvature perturbation 
on uniform density hypersurfaces. When the curvaton 
decays into radiation, the total adiabatic perturbation is 



C = r a ( a + (1 - r a )( inl 

— Cinf + -^S a 



(3) 



and "inf" denote the curvaton and inflaton 
perturbations, respectively, and 



3» g 
l-O, 



(4) 



gives the relative contribution of the curvaton to the to- 
tal curvature perturbation at the time of decay. S a is 
the non-adiabatic part of curvaton perturbations and is 
defined as S a = 3(C CT — Cm/)- From this expression it is 
evident that the contribution of the curvaton grows as it 



oscillates in the post-inflationary universe before decay: 
an oscillating scalar evolves as p a oc a~ 3 , and radiation 
as p 1 oc a~ 4 , giving Pa/p^ ~ ^er oc a. 

Since the adiabatic inflaton perturbations and the 
isocurvature curvaton fluctuations are uncorrelated, the 
total power spectrum of primordial curvature perturba- 
tions is given by 



P(k) 



fc 3 ICI 2 

2tt 2 

Pinf(fe) 



-PsAk)- 



(5) 



One can therefore model the total power spectrum of 
the mixed curvaton-inflaton perturbations to lowest or- 
der in slow-roll, as a combination of two power laws 



P(k) = A s 



(1-/) 



Tlinf — 1 



which is of the form Eq. (TTJ) with 

Pni{k ) 



A. 



f 



1-/ ' 

9P mf (k )+r 2 .P Sa (k Q )- 



(6) 

(7) 
(8) 



The function / has been introduced to control the rela- 
tive contributions of the two spectra, allowing the com- 
posite spectrum to range from pure curvaton (/ = 1) to 
pure inflaton (/ = 0). In the absence of additional dis- 
criminating observables like isocurvature modes or NG, 
it is generally not possible to separately constrain the 
curvaton and inflaton contributions; the resulting degen- 
eracy is not easily reduced [TTJ. We emphasize that the 
assumption of a pure power law for both the inflaton 
and curvaton components is restrictive: in principle, ei- 
ther spectrum could also have intrinsic scale dependence, 
which would add additional parameters to the fit. 

The effective spectral index and running of the full 
bimodal spectrum Eq. (pi) can be obtained by Taylor 
expanding the spectrum in Ink. From the definitions 



dlnfc 


(9) 


d 2 lnP(fc) 
dlnfc 2 _ "' 


(10) 


we obtain the relations, 




a 


(11) 


n s 1 ^inf 1 -. n > 


a = (l-/)/(An) 2 , 


(12) 



where An = rii n f — n a . One notable feature of the spec- 
trum is that the effective running is positive, a > 0, 
a simple consequence of the fact that the component 
with lower spectral index will always dominate on large 
enough scale, and the component with higher spectral 




FIG. 1: Depiction of a as a function of An and /. Running 
is largest when both An and / are large. 



index will always dominate on small scales. This is in 
tension with current cosmic microwave background mea- 
surements which indicate a preference for a < [12H14J . 
This fact can be used to constrain the bimodal spectrum, 
as a model, relative to the power law spectrum of the pure 
curvaton or inflaton. And, while measurements of n s and 
a are not sufficient (by themselves) to determine the un- 
derlying parameters n in f, n CT , and /, Eqs. (11) and (12) 



can be used with precision measurements of n s and a to 
establish useful constraint relations between the underly- 
ing parameters. Additional observables, like isocurvature 
modes or NG, can be used to further improve the recon- 
struction of the bimodal spectrum if they are observed 
in the future. We consider this possibility as well. 



III. CONSTRAINTS FROM CMB 
MEASUREMENTS 

In this section we obtain constraints on the bimodal 
power spectrum Eq. (|6J). It is assumed that ri[ n { and 
n CT share the same prior range so that / need only vary 
between and 0.5. This restriction gives the parameter- 
ization 



Pn(k) = A s 



(1-/) 



koj \k 0/ 

(13) 
where "dom" labels the dominant and "sub" the sub- 
dominant spectrum, measured at a pivot scale of ko = 
0.017 Mpc" 1 h. We use the 7-year WMAP dataset 
(WMAP7) in combination with data from the South Pole 
Telescope (SPT) [T2] Q3] to perform a Bayesian analysis 
of the power spectrum in terms of the underlying parame- 
ters A s , f, n dom , and n su b- The base cosmological param- 
eters are chosen to be the baryon and cold dark matter 
densities, il^h 2 and Vt c h 2 , the ratio of the sound horizon 



to the angular diameter distance at decoupling, 9 S , and 
the optical depth to reionization, r. We also vary the 
tensor/scalar ration, r. Use of the SPT data introduces 
three more foreground parameters: the Poisson point 
source power from randomly distributed galaxies, Df°, 
the clustered point source power, Df , and the secondary 
emission from the Sunyaev-Zeldovich effect from clusters, 
Asz- We also obtain constraints for a simulated Planck- 
precision data set that we describe below. We perform 
the parameter estimation using Markov Chain Monte 
Carlo (MCMC) with CosmoMC [15], marginalizing over 
the foreground parameters. Chain convergence was de- 
termined with the Gellman-Rubin statistic, R— 1 < 0.1. 
While ground-based observatories like the Atacama 
Cosmology Telescope (ACT) and SPT have provided in- 
creasingly precise measurements of a, ESA's Planck Sur- 
veyor satellite promises the best constraints yet on both 
n s and a. We therefore obtain constraint projections 
on the bimodal spectrum based on a simulated Planck- 
precision data set in addition to WMAP7+SPT. We per- 
form MCMC on the same set of base cosmological pa- 
rameters as above: flbh 2 , fl c h 2 , 9 S , r, with spectral pa- 
rameters A s , ndom, "-sub, and /, using simulated tem- 
perature and polarization data out to I = 2000. We 
conduct the analysis on a fiducial model with a power 
law spectrum and r = 0. The base cosmological param- 
eters (ilbh 2 , £l c h 2 , Sl t, A s , and n s ) are taken equal 
to the WMAP7+SPT best-fit values. Data is simulated 
for three channels with frequencies (100 GHz, 143 GHz, 
217 GHz) and noise levels per Gaussian beam (cj ix ) 2 = 

(46.25 fiK 2 , 36 /xK 2 , 17.6 fiK 2 ) (with <t£ x = %/2crJ ix for 
polarization spectra). The FWHM of the three channels 
are #f w hm = (9.5', 7.1', 5.0') and we assume a sky cover- 
age fraction of 0.65 [15] . 

Although the MCMC furnishes constraints directly on 
the base parameters of Eq. (13), we first examine es- 
timates for the derived parameter n s . We compare con- 
straints in the n s -r plane obtained for a power law ansatz, 
P(k) ex fc" s_1 , with those obtained for the bimodal spec- 
trum Eq. (13). The bimodal spectrum can be written 



as a single spectrum with an effective spectral index and 
running, obtained from the base parameters via Eq. ( 11 ). 



The bimodal spectrum therefore effectively has one more 
free parameter than the power law, and in general, the 
addition of a free parameter can improve a model's fit to 
data simply on account of the reduced number of degrees 
of freedom. If the additional parameter has correlations 
with the base parameters, degeneracies might open up 
and constraints on the correlated parameters might de- 
grade. This is what happens in the n s -r plane when run- 
ning is added as a free parameter. We present results in 
Figures [2] and [3] for the three data sets: (a) WMAP7, (b) 
WMAP7+SPT, and (c) Planck. Interestingly, none of 
the constraints obtained here suffer from any such degra- 
dation, with the contours in each case virtually overlap- 
ping. This is because the additional effective parameter 
of the bimodal spectrum is positive running, which is not 
well-tolerated by current CMB data: -0.060 < a < -0.008 
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FIG. 2: Marginalized constraints in the effective n B -r plane for the bimodal spectrum for (a) WMAP7 and (b) WMAP7+SPT. 
For each data set we also include results from a power law fit (dashed lines). The two sets of contours overlap in each case. 
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FIG. 3: Marginalized constraints in the effective n s -r plane 
for the bimodal spectrum for a simulated Planck-precision 
data set, showing a power-law fit (dashed lines). 



(WMAP) Q2] and -0.037 < a < -0.011 (WMAP+SPT) 
[2] at 68% CL, both assuming a zero tensor amplitude. 
Constraints are weaker when r is included as a joint pa- 
rameter in the fit [T7] . For the bimodal model, the points 
of appreciable likelihood - those that contribute to the 
contour - are those for which a ~ 0. At the pivot scale 
(fco = 0.017 /iMpc - ), positive running forces a red spec- 
tral index. A red spectral index with positive running 
will fit the data worse than a red spectral index with no 
running, because the positive running results in an ex- 
cess of power at the low multipoles and a blue spectrum 
at high multipoles. These features are not in agreement 
with CMB data IT51. 



This fact is also revealed by examining the overall 
model fits: the bimodal spectrum gives a best-fit log- 
likelihood — 21n£ = xiff — 7514.25, while the power 
law gives xlff = 7514.4 for WMAP7+SPT, with simi- 
lar close agreement for the other data sets. The reduced 
Xeft, given by xlff/v, is the relevant statistic for model 
comparison. Here v denotes the number of degrees of 
freedom, defined as the number of data points minus the 
number of free parameters. With two fewer degrees of 
freedom (i.e. two additional parameters), the bimodal 
spectrum gives a larger reduced Xeff an d is consequently 
a worse fit to the data. Power law spectra, whether due to 
a pure inflaton or curvaton, are therefore preferred over 
the bimodal spectrum arising from mixed perturbations. 

Goodness-of-fit aside, current data does not go so far 
as to rule out the bimodal power spectrum. It is there- 
fore of interest to examine how well CMB measurements 
constrain the parameters describing the bimodal power 
spectrum, since these quantities are directly related to 
the underlying inflationary degrees of freedom. As men- 
tioned, the base parameters n^ omi n su b, and / cannot be 
separately determined from the overall power spectrum, 
but quality measurements of n s and a can establish con- 
straint relations among them. From Eq. (12), a quality 



measurement of running will furnish a relation connect- 
ing An = ridom — «sub and the fraction /. We present 
results in the An-/ planes in Figures [4] (a) , (b), and (c). 
The distributions are shaped by the quality of the con- 
straint on positive running: the inclusion of SPT data 
improves the constraint relative to WMAP7 alone, 1 with 



1 It is notable that the improvement in constraints in Figure |4J(b) 
relative to (a) is less striking than would be naively expected 
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FIG. 4: Marginalized constraints in the An-/ plane for (a) WMAP7, (b) WMAP7+SPT, and (c) projection for a Planck- 
precision CMB data set. 



Planck offering the highest precision. With Planck, it be- 
comes possible to obtain constraints on the fraction /: for 
very different spectral indices ndom and n S ub, the fraction 
is constrained to be very small. Conversely, for small /, 
the spectrum determined by ridom completely dominates 
and n su b is unconstrained. For / = 0.5, constraints are 
symmetric about An = 0. Meanwhile, for intermediate 
/, An < means that the dominant spectrum is red- 
der at the pivot, with the subdominant spectrum coming 
to dominate on smaller scales, leading to a milder over- 
all positive running. For intermediate / and An > 0, 
the dominant spectrum is bluer at the pivot than the 
subdominant spectrum, which might therefore come to 
dominate on large scales, again leading to a mild overall 
positive running. These two cases should not, however, 
be equally favored by the data owing to the fact that 
SPT places strong constraints on the high-£ spectrum, 
while the low-Fs remain poorly constrained due to cos- 
mic variance. There appears to be a slight indication of 
this asymmetry in our results, most pronounced in Figure 
H (c) where constraints on An are not symmetric about 
zero, and tighter on An < 0. 

While current data provide no evidence for mixed per- 
turbations, and while future data sets will begin to re- 
solve the bimodal spectrum, the inflaton/curvaton de- 
generacy remains unbroken. It is well known that, in 
general, additional observables beyond the adiabatic and 
tensor power spectra on CMB scales are needed to re- 
solve the degeneracy. It is shown in [TT] that a detec- 
tion of isocurvature perturbations would not only break 
the degeneracy but would give the correct number of ob- 
servables to successfully resolve the parameters of the 
underlying model. In particular, a measurement of the 
correlation angle between the adiabatic and isocurvature 



components, 



given the rather substantial improvement in constraints on a 
from SPT. The effect is one of priors: constraining the bimodal 
spectrum is equivalent to fitting a power-law with positive run- 
ning; the prior on a corresponds to a region of low relative like- 
lihood where a is not sharply constrained 



cosA 



Pcs(k ) 



^/p c (k )p s (k y 



(14) 



will provide sufficient information to determine the frac- 
tion /. The curvaton can generate an isocurvature per- 
turbation in a variety of ways, and the degree of cor- 
relation between the isocurvature and adiabatic modes 
depends on the process. For example, a cold dark mat- 
ter isocurvature mode can be set up if the dark matter 
is a direct product of curvaton decay. In the pure cur- 
vaton limit, this mode is anticorrelated (/3 = —1) with 
the adiabatic mode. The correlation angle can be written 



I3 = - 



1 



Vl + A- 



(15) 



where A parameterizes the contribution of the curvaton 
to the overall spectrum, Pq = (l + A)P; n f. The parameter 
A is related to the fraction / through 



/ 



A + l' 

2 



P 



(16) 



The degree of correlation varies with the fractional con- 
tribution of the curvaton to the overall spectrum: for 
anticorrelated modes, / = 1 as expected. While current 
CMB data shows no evidence for an isocurvature mode, 
it is still possible that isocurvature perturbations repre- 
sent a small fraction of the overall density perturbation. 
If an isocurvature mode is detected by Planck, the corre- 
lation angle is expected to be measured with a precision 
of 5/3 ~ 0.04 at 68% CL. This leads to a constraint of 
8f = 2/3(5/3 on the curvaton fraction, which, from Figure 
3 (c) for / > 0.1, gives |An| w 0.1 at 68% CL. 

Another possibility of resolving the innation/curvaton 
degeneracy is to detect a significant level of NG. If 
we parametrize the NG by the dimensionless parameter 
/nl ~ (C 3 )/(C 2 ) 2 an d assume that NG is sourced only by 



the curvaton, we find 



/] 



NL 



(17) 



The Planck experiment is sensitive to values of /nl as 
small as 0(5). In the case in which no NG is measured, 
we conclude that r a > / 2 /5. A possible detection of a 
contribution to the power spectrum from a second con- 
tribution at the 10% level would then require r a > 0.002. 



IV. CONCLUSIONS 

In this paper, we consider CMB constraints on a mixed 
inflaton/curvaton power spectrum, which we model as a 
superposition of two power laws with independent spec- 
tral indices. Locally, such a bimodal spectrum can be 
well approximated by scale-dependent power spectrum 
with positive running a, since the redder (smallest n) 
spectrum will dominate at large scale, and the the bluer 
(largest n) spectrum will dominate at small scale (Fig. 
[l]). This is significant in the context of current data, 
which favor either a pure power law or negative running 
[L2TfT4*] . Consistent with this expectation, we find that 
adding new parameters describing an effective positive 
running does nothing to improve the quality of the fit. 
In particular, constraints on the effective spectral index 
n s vs. tensor/scalar ratio, r, are almost completely un- 
affected by the addition of the new parameters (Figs. [2] 
pi), consistent with a single power-law spectrum being fa- 
vored by the data. In principle, either the inflaton or the 



curvaton spectrum could have an intrinsic scale depen- 
dence, which would complicate the analysis relative to 
simple assumption of two pure power laws assumed here. 
This would further reduce the discriminating power of 
the data. 

The parameters /, specifying the relative amplitudes 
of the power spectra, and the difference in spectral 
indices An, are poorly constrained (Fig. H). This is 
true both for the WMAP7 data considered alone and 
in combination with South Pole Telescope. In the case 
of the simulated Planck data, we find that Planck will 
be able to significantly constrain the fraction / in a 
subdominant spectrum as long as the spectral indices 
differ by a sufficient amount, |An| > 0.2. In this case, 
a contribution to the power spectrum from a second 
power law can be ruled out at better than the 10% level 
at 95% confidence by Planck (Fig. g). For |An| < 0.2, 
constraints from Planck will be poor, which is to be 
expected, since in the limit that both power spectra 
are identical (|An| -c 1) the spectrum is effectively no 
longer bimodal, and no constraint on / is possible. 
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